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Abstract 

Field theories which violate the null energy condition (NEC) are of interest 
for the solution of the cosmological singularity problem and for models of cos- 
mological dark energy with the equation of state parameter w < — 1. We discuss 
the consistency of two recently proposed models that violate the NEC. The ghost 
condensate model requires higher-order derivative terms in the action. It leads 
to a heavy ghost field and unbounded energy. We estimate the rates of particles 
decay and discuss possible mass limitations to protect stability of matter in the 
ghost condensate model. The nonlocal stringy model that arises from a cubic 
string field theory and exhibits a phantom behavior also leads to unbounded en- 
ergy. In this case the spectrum of energy is continuous and there are no particle 
like excitations. This model admits a natural UV completion since it comes from 
superstring theory. 
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1 Introduction 



There are general restrictions to the energy-momentum tensor T Mi , of a physical system. 
Such restrictions are referred to as energy conditions. They play an important role in 
general relativity, in particular in considerations of the black holes and cosmological 
singularities [TJ [2] (see also [3]- [7] and refs. therein). A weak form of the energy 
condition states that T tlu n fJ- n u > for any null vector n M ; it is called the null energy 
condition (NEC). 

The NEC is directly related with a restriction on the dark energy equation of 
state parameter w = p/p. The energy-momentum tensor of dark energy is T^ u = 
diag(p,p,p,p) with positive energy density p and negative pressure p [HE]- The con- 
dition w < — 1 implies violation of the NEC. Since experimental data allows a 
possibility of w < — 1 the study of such models attracts a lot of attention. Models of 
bouncing from cosmological singularity also require violation of the NEC [2l [TTj . 

There are general results that coupled scalar-gravity models which violate the NEC 
are unstable [3l SI El d EE2l IS] but recently there have been proposed the ghost con- 
densate model [13l [HI [151 IS] an d the nonlocal stringy model [T71 [19j [18j [20] which 
apparently lead to consistent effective theories. Both these models include higher-order 
derivative terms. In this note we analyze these two models using the Ostrogradski 
method. 

The ghost condensate model requires higher-order derivative terms in the action 
hence it leads to a massive ghost field and unbounded energy and as a result it could 
violate the stability of matter. We estimate the rates of particles decay and discuss 
possible mass limitations to protect stability of matter in the ghost condensate model. 

The nonlocal stringy model that exhibits a phantom behavior also leads to un- 
bounded energy but in this case the spectrum of energy is continuous and there are no 
particle like excitations. This model admits a UV completion since it comes from the 
superstring theory. 

In more details, the model is based on SFT formulation of a fermionic NSR string 
with the GSO— sector [21J. In this model a scalar field is the open string tachyon 
|22j . which describes according to the Sen conjecture [23J a dynamical transition of a 
non-BPS D-brane to a stable vacuum (see [23] for review). This stable string theory 
is supposed to be described by a VSFT (Vacuum String Field Theory). The model 
that we are going to present is an approximation to VSFT which is stable and has no 
particle excitations corresponding to the open string. 

A vector-scalar model violating NEC with excitations modes stable in some region 
has been proposed in [25] . There is a tensor-scalar dark energy model admitting phan- 
tom behavior at small redshifts|26j. String-inspired and braneworld dark energy models 
are also the subject of intensive study in the last years (see for example [27J [281 I2H] 
and refs. therein). 

The paper is organized as follows. In Section 2 we discuss the ghost condensate 
model. We show that the model actually includes two non-relativistic scalar fields. One 
field is a massive ghost while the other field has a positive energy. We evaluate the 

1 A direct search strategy to test the inequality w < — 1 has been proposed [ID] . 



2 



decay rate of an ordinary particle in the ghost condensate model and obtain a rather 
stringent limit to the mass of the particle if the lifetime is greater then the Hubble 
time. In Section 3 we consider the nonlocal stringy model. We used the Weierstrass 
product representation and the Ostrogradski method to display the spectrum of the 
model. We obtain that also for this model the energy is unbounded from below but in 
this case there is no particle like excitations. 



2 Ghost Condensate Model 
2.1 Setup 



The ghost condensate model in Minkowski space (signature is (— , +, +, +)) deals with 
the Lagrangian [131 CGQ E51 EH] 

C = M 4 P(X) + M 2 S!(X)(D0) 2 + M 2 S 2 {X)d^d u ct)d^d u (t) (1) 

Here is a scalar field, 

X = -d^(j> = 2 - (V0) 2 (2) 

and M an arbitrary mass scale. Ghosts (phantoms) could come from the term P(X). 
There are higher-derivative terms in the Lagrangian proportional to Si and S2 which 
also lead to ghosts as it will be discussed below. 
If 

P(X)= l -{X-lf = --^ + ... (3) 

then the kinetic term has the "wrong" sign and one has ghosts. This corresponds to 
the "wrong" vacuum state solution of the field equations 

0o = (4) 

It was proposed in [T3] that to cure the theory, by analogy with the Higgs mecha- 
nism, one can consider the configurations with P'{X) = as a candidate ground state. 
There is a solution of the field equations of the form 

0o = t (5) 

which satisfies P'(X) = 0. 

If we consider small fluctuations n(t, x) about this solution 

= t + 7r(t,x) (6) 

then the Lagrangian ([3]) for quadratic fluctuations leads to the kinetic term with the 
"proper" sign 

C = -M 4 7T 2 (7) 
2 v ; 
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There is no the spatial kinetic term (Vir) 2 in the quadratic Lagrangian for tt. The 
higher order terms proportional to Si and S2 in ([1]) are added to P to get a spatial 
kinetic term. The following quadratic Lagrangian is considered in [131 EE HH IB] 



£ = I M 4 7r 2 -i 7 2 (V 2 7r) 2 
2 2 1 y ' 



where 7 2 is a constant, 



7 2 = -2M 2 (S 1 (l) + S 2 (l)) (9) 

We note, however, that the original Lagrangian ([T]) leads not to ([H]) but to the 
following quadratic Lagrangian 

C = Vtt 2 - i 7 2 (D7r) 2 (10) 



where 



□ TT = -ft + V 2 7T. 



The Lagrangian ffTUj) includes the higher derivative term which leads to ghosts for all 
momenta. In fact it describes not a single but two scalar fields. To demonstrate it let 
us consider theories with higher derivatives. 



2.2 Ostrogradski Method for Equations of Higher Order 

Let us consider the Lagrangian 

C = <pn(a + /3n)(p (11) 

where a and (3 are real parameters. By using the Ostrogradski method [3Q1 EI] we 
introduce the fields 

$ = (a + pn)<p, <f> = D<p (12) 
Then, modulo the surface terms, the Lagrangian (TTTT) can be written as 

£ = -ijjUif)- P(f)( a + f3n)(J) (13) 
a a 

We see that independently of the signs of a and (3 one of fields and if) is a ghost. 
For instance if a > and (3 < then we have the massless ordinary field if) and the 
massive ghost <f>. 
The Lagrangian 

£ = <pn 2 ip (14) 

is more complicated for studies. The limit a — > is singular and to compute it one has 
first to make a canonical transformation [32, 33J. 
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2.3 Ghosts in the Ghost Condensate Model 

Now let us discuss the Lagrangian (jTOj) 

2 2 ' v ; 

Equation of motion is 

(□ 2 + f/d^TT = (15) 

where /z 2 = M 4 /7 2 or 

(5f - 2<9 2 V 2 + V 4 + f/d^n = (16) 
For the spatial Fourier transform of the field n(t, k) the equation (JT6]) reads 

(df + (2k 2 + fi 2 )d 2 + k 4 )n = (17) 

The characteristic equation 

u; 4 - (2/c 2 + /i 2 )cg> 2 + A; 4 = (18) 

has only real roots 

< 2 W = ^ 2 [l + 2^±Wl+4^] (19) 

and we obtain 

df + (2k 2 + ^ 2 )d 2 + k* = (d 2 t + uj 2 )(d 2 + u 2 2 ) (20) 

Denote 

I ~fc2 

a 2 = cj 2 -cj 2 = y u 2 Wl + 4— >0 (21) 

By using the Ostrogradski method we introduce the fields 

^ = (d 2 t+ ul)-K, <P = (d 2 t + u 2 2 )-n (22) 
a a 



and obtain that the Lagrangian ([10 

1 



C = -- 7 2 vr[9 4 + (2A; 2 + /i 2 )«9 2 + fc 4 ]vr (23) 

modulo the surface terms can be represented as 

£ = \i 2 <P(d 2 + ujf)^ - i 7 V(«9 2 + u 2 > (24) 

We obtain two non-relativistic fields <fi and ip. If r y 2 > then the field <fi leads to ghosts 
for any momenta k. Note that there is a mass gap for ghosts since u>i(k) 2 > fi 2 . For 
large k one gets a relativistic dispersion law 0J\(k) ~ a/A; 2 + /i 2 . The field ?/> does not 
have a mass gap. 
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2.4 Decay of Particles 



Ghosts (phantoms) could have negative energy. Therefore ordinary particles can decay 
into heavier particles plus phantoms. This is discussed in J3JH]. 

We want to estimate the lifetime of the ordinary particle under this decay. We do not 
have a direct interaction describing this decay. What we have for sure is an interaction 
of the phantom with gravity and interaction of gravity with the usual particles. A 
simple diagram describing a decay of an ordinary particle ijj into usual particles ip and 
X and two phantoms is presented in Fig.l. There h is a graviton. 



Figure V. Decay of ip particle into two phantoms and two other particles ip and \ 

We suppose that the interactions "gravity-phantom" and "gravity-ordinary fields" 
come from the mass terms in the action 



Here is a phantom and tp and \ are ordinary fields. There is a cubic interaction 
of the ordinary fields but they couple with the phantom only through gravity. We 
assume that the mass M x of the field x is °f the order of the large phantom mass M 
and the mass of the field ip is much smaller. From this action we get the following 
interaction of gravity with matter fields 




v. 



S = M 2 J d A x^gR +-J d 4 xy^g[{d(f)) 2 + M 2 2 
-(cty) 2 - m 2 ^ 2 - (di;) 2 - M 2 X 2 + \^ 2 x\ 



(25) 




(26) 



(27) 



where h is a scalar field symbolically describing graviton. Therefore, the coupling 
constants of gravity with the and ip fields are 



The decay rate of a particle to phantoms through the channel (see Fig.l) 

^->^ + X + 0i + 02 (29) 

is 



J_ f A d 3 P(f>l d 3 p 4>2 d 3 Pi , d 3 p x 2 
m^J (27r) 3 2E 4>1 (27r) 3 2E <h (2n) 3 2E lP (2n) 3 2E x lJ 1 1 ] 

(27r) 4 5 {4) (p^ n - P<h - P(/>2 -p^p- p x ) , 
where the matrix element M. is 

1 1 

M = A^A^Aw^ — — — — - (31) 
p 2 q 2 + 

Here A is a momentum cutoff, p is the momenta of the internal graviton line and q of 
the internal ip line. We assume that the dimensional coupling constant A x ^2 is of the 
order 

Using (EHD and flZZ} we get 

m^M 2 1 1 
^ = ~Mf7q 2 + ml- (32) 



The integral is convergent even without cutoff. To estimate it we can take 
We estimate the decay rate as 



M ~ M p ~ A 



(33) 

The timescale for decay is = 1/r. A model will be phenomenologically viable if the 
lifetime is greater than the Hubble time Hq 1 , i.e. 

5 

T ~ M p <H ~ l(T 60 M p (34) 

p 

Using M p ~ 10 19 GeV we get 

< 10 7 GeV (35) 

A stringent limit to the mass one can obtain from the investigation of the quartic 
interaction of the ordinary particles of the form \ip 3 x with the coupling constant A of 
the order one. Then instead of ( 1331) and (1351) we get 

m 3 ,, 

r - Jfs M P ( 36 ) 
p 

and an uncomfortable result 

< 10 _1 GeV. (37) 
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3 Nonlocal Stringy Model 

3.1 Example of Field Equations of the Exponential Type 

Consider first the following action 

£ = 0( e D - 1)0 (38) 

The field equations are nonlocal since one can write them as integral equations. It 
serves as an example for the nonlocal stringy Lagrangian which will be considered 
below. We don't consider here the rigorous mathematical approach to such equations, 
see recent papers [34] . By using the Weierstrass product 

oo 2 

e *_l = e */2 iZ JJ(i + £ ? ) (39) 

where 



uj 2 . = 4ir 2 j< 

we write 



n 2 

e n_ 1 = e n/2 D TT (1 + ^ ) (40) 
Therefore [32] the Lagrangian fl38l) has the same spectrum as 

oo 

i=i 

Here 

£j = l/^F'{-^) (42) 

where 

, 9s . , , e~ z l 2 sinhz/2 , , 

F(z 2 ) = (e z -l) = 2 '— (43) 

z z 

More general Lagrangian 

C = <f)(e an - \)(f> (44) 

appears in p-adic string theory (see [33 EE] and refs therein). Here a and A are 
constants. The spectrum of the theory can be read out of the Weierstrass product 

e aD - \ = X^e an /\aD - log A) TT(1 + ^£^^1!) (45) 
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3.2 Nonlocal Tachyon Field 

An effective Lagrangian for the tachyon field coming from GSO— sector of the fermionic 
NSR string which is described by cubic string field theory [21] has the following form 

£ = l0 e -3 n (n + m 2 )0-^0 4 (46) 
There is the "wrong" vacuum state solution of the field equations 

0o = (47) 
which leads to the tachyon equation for fluctuations about this vacuum: 

(□ + m 2 )0 = O (48) 
There exists also another vacuum solution 

0o = m/V\ (49) 
One obtains the following Lagrangian for the small fluctuations about this solution: 

2 

Tf) 

C = — 0F(-D)0 (50) 

where 

F(*) = (-£ 2 * + 1) e** - 3, (51) 

and we set £ 2 = 1/m 2 . 

The characteristic equation 

F{z) = (52) 

has an infinite number of complex conjugated pairs of complex roots z — k? (see J2U] 
for the consideration of these roots in the case corresponding to the action (l5"0l with 
( |5T|) and [37] for the case of arbitrary constants). The roots n 2 are expressed by means 
of the Lambert function W(z) satisfying the equation W(z)e w( - Z,> = z. 
We represent (!5Tl) in terms of the Weierstrass product 

F{z) = -2e^ 2 -^ I^ 1 - - 4) ex PK \ + 4)] (53) 

Now we can write an appropriate quadratic Lagrangian (up to total derivatives) as 

^ = E e ^( D - K ?)( D - K f)^ ( 54 ) 



Equations of motion with complex distinct frequencies were studied in [32]. Con- 
sider the following equation: 

(□ - k 2 ) (□ - k* 2 ) = (55) 
9 



where k = v + ia. The corresponding Hamiltonian is a linear superposition of pairs of 
complex conjugated oscillators of the form: 

H=\ + ^(kfQl(k) + P 2 {k) + u(kr 2 Q 2 2 (k)) (56) 

k 

where 

P\ — P'2, Q* = Q2 

and 

Uj(k) 2 = K 2 + k 2 

The spectrum of the Hamiltonian has the form 

nv(k)+pa(k) (57) 

where we put oj(k) = u(k) + ia(k). Here n — 0, ±1, ... while p is continuous and ranges 
from — 00 to +00. The energy spectrum is indefinite and continuous. Therefore no 
particle attributes can be ascribed to these modes [52] . 

4 Conclusions 

We have considered two recently proposed models which violate the NEC. In both 
models there are higher derivative terms and we have used the Ostrogradski method to 
study them. A straightforward application of this method to linear theories indicates 
that energies are unbounded. Note that the problems with unbounded energy consid- 
ered in this paper are mathematically similar to the eigenvalue problem for the non 
positively defined hyperbolic Klein-Gordon equation on Lorentzian manifolds which is 
solved in [38] . 

One can expect that an incorporation of nonlinear terms could drastically change 
the situation. In particular, in nonlinear theories could exist islands of stability [H 1391 

annul. 

In the ghost condensate model with higher-order derivative terms it would be in- 
teresting to find a UV completion of the theory and study islands of stability. 

The nonlocal stringy model has been proposed to describe a decay of D-brane and 
though it also leads to unbounded energy in this case the spectrum of the energy 
is continuous and there are no particle like excitations. This model admits a UV 
completion since it comes form superstring theory but it would be very interesting to 
find a direct mechanism of compensation of unbounded continuous spectrum. 
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